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Abstract 

Divergences that arise in the quantization of scalar quantum field 
models by means of a lattice-space functional integration may be at- 
tributed to a single integration variable, and this fact is demonstrated 
by showing that if the integrand for that single integration variable is 
appropriately changed, then a perturbation expansion becomes order- 
by-order finite and divergence free. The paper concludes with a brief 
review of a current proposal of how an auxiliary, nonclassical poten- 
tial added to the lattice action of a relativistic scalar field quantization 
may automatically render an analogous change of the integrand, and 
thus may lead, as well, to nontrivial and divergence-free results. 

1 Introduction and Overview 

Traditional formulations of quantum field theory often encounter one or an- 
other divergence in the course of calculation. For a self-interacting, relativis- 
tic scalar field, for example, the divergences are often classified on the basis 
of those that are encountered in the course of a perturbation analysis. If we 
let n denote the number of space-time dimensions, then a typical classical 
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action functional is given by 

I = J {|[^0(a;)] 2 - \m 2 (f)(xf - g(j)(x) p } dhc 

appropriate to a p th power interaction term. This is the action for the so- 
called tp v n model, and we shall confine attention to such models. Here, there 
is a sum over the indices \x from to s — n — 1, and the metric signature 
is [+1, — 1, — 1, . . . , — 1]. The kind of divergences that arise in quantization 
depend strongly on the space-time dimension n. For example, for n — 2, the 
only divergences that arise have to do with normal ordering which amounts 
to a (re) definition of the local product for operators. For n > 3 the kind of 
divergences includes those covered by normal ordering and generally addi- 
tional divergences as well. In particular, whenever p < n/(n — 2) it follows 
that normal ordering cures all divergences on a term by term basis within a 
perturbation expansion. When p < 2n/(n — 2) we are in the so-called super- 
renormalizable regime, in which case infinite renormalization counterterms 
are generally required, such as a divergent mass renormalization counterterm 
for the ip\ model. When p = 2n/{n — 2) we are in the strictly renormaliz- 
able case where there is generally both coupling constant and wave function 
renormalization. Finally, when p > 2n/(n — 2) we are in the nonrenormal- 
izable situation and the number of divergent interaction counterterms that 
need to be added grows without limit. This brief survey indicates the normal 
situation viewed from the point of view of perturbation theory. For further 
details see [0. 

One may also formulate things differently. If we choose a Euclidean lattice 
space functional integral formulation, then on the basis of the renormalization 
group, for example, and outside of perturbation theory, one can show that 
the nonrenormalizable cases, i.e., where p > 2n/(n — 2) actually pass to a free 
theory, or at least a generalized free theory, in the continuum limit. Thus 
for the very important case where p — 4, it follows that for n > 5 the theory 
becomes trivial in the continuum limit ||. Although the full proof has not yet 
been established there is every indication that the same conclusion actually 
holds in the case p = 4, n = 4. Thus we observe that different calculational 
schemes may very well lead to distinct answers relative to the question of 
divergences. 

Let us recall the appearance of the Euclidean-space lattice formulation 
for the cpP model. We let k = (fc , h, ■ ■ ■ , K) , kj e Z = {0, ±1, ±2, . . .}, 
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denote a lattice site, k* denote any one of the nearest neighbors (in the 
positive direction) to the point k, and a denote the lattice spacing. Then 
the generating function for the Euclidean-space correlation functions on the 
lattice is given by 

S(h) =C N J exp{ m k <p k a n - \Y{<j) k * - (j) k ) 2 a n - 2 

-\ml^\a n - g^la n }Iid<p k . 

This integral is to run over N = (2L + l) ra lattice sites of a (hyper)cubic 
lattice, and the continuum limit is one in which a — > as well as L — > oo 
in such a way, ultimately, that the volume of the space-time integration 
region becomes all of R n . The constant CV is fixed by the requirement that 
S(0) = 1. Here we have introduced unknown parameters which are to be 
determined, if possible, so that a consistent and nontrivial (non- Gaussian) 
theory emerges in the continuum limit. It is this formulation which, even 
allowing for arbitrary choices of the parameters Y > 0, m^, and go > 0, 
nevertheless results in a trivial, Gaussian behavior in the continuum limit, 
i.e., when a — > 0, for n > 5 (n = 4). These results emerge from an application 
of the renormalization group and/or the random- walk representation 0]. 

In the present paper we wish to sketch an argument which leads to a 
qualitative appreciation of various divergences and which at the same time 
isolates a single parameter which, in a manner of speaking, can be said to be 
'responsible' for all the divergences. This identification is significant because 
in work reported elsewhere — and very briefly summarized in the final sec- 
tion of the present paper — we have shown how one may effectively eliminate 
an analogous parameter in lattice formulations of more physical models. This 
procedure, when carried out to its potential — a task, unfortunately, that has 
yet not been fully accomplished — has the capability of offering a quantum 
field theory the perturbation theory of which is order-by- order divergence 
free! 

2 Hyper-extreme Spherical Coordinates 

The lattice-space functional integral involves an integration over iV fields on 
the N distinct lattice sites. Let us reexpress that integral in new coordinates 
which we shall refer to as hyper- extreme spherical coordinates. To that end, 
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consider the measure 

IT d<p k = IT d<j> h 5(k 2 - T.cpl) dn 2 = 5(k 2 - rf/t 2 n 

= k n ~ 2 dn 2 5(1 - £?? 2 ) n d Vk = 2k n ~ 1 dn 5(1 - Y.rg) U drj k . 

In the last line we have introduced <p k = Ki] k for all sites k in the lattice. 
Observe that k represents the radius of all the N fields while {r/ k } denotes 
the direction field of the N field amplitudes {<f) k }. The direction field satisfies 
Hrjl = 1, and thus this field lies on (the surface of) the (N — l)-sphere, 
S*^ -1 . Let us introduce these variables into the lattice regularized functional 
integral. It follows (with Cjv absorbing the factor 2) that 

S(h) =C N J exp{Km k r] k a n - i/t 2 FS( % . - r] k ) 2 a n - 2 - \m 2 K 2s Lr) 2 a n 

-g K p Zr] p k a n } k n ~ 1 dn 5(1 - S^ 2 ) n drj k . 

Now observe that there are two qualitatively different kinds of integration 
variables. On the one hand, there are the r] k variables constrained so that 
the sum of their squares lies on the unit sphere. Such variables are quite 
'tame' in their behavior, and as we shall see they are not responsible for the 
divergences even though there are a great many of them. On the other hand, 
there is the single variable k which enters into the integrand in one factor 
with a large power, namely the term /t^" 1 . If we look at the behavior of the 
k integral we see that it is largely controlled by that term with a high power 
(N — 1). A steepest descent integration for k is suggested and leads to the 
fact that k oc y/N whenever p and n satisfy p < 2n/(n — 2) (see below). Not 
only does the steepest descent approach lead to a large value for k, it also 
leads to a relatively very narrow width of order one. Of course, a change of 
variable given by rescaling k could put the new support of order one, but 
then it would be extremely sharply concentrated. Essentially, it is the fact 
that the support of k is extremely sharply concentrated that ultimately leads 
to a divergent perturbation expansion. Let us illustrate and thereby clarify 
this remark. 

Consider first the case of a free field in n dimensions; specifically we 
consider 

S(h) = C N J exp{Km k r] k a n -\K 2 ^(n k * -n k ) 2 a n ~ 2 -\m 2 K 2 a n } k n ~ 1 dndairj) 
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where we have set Y — 1 and mo = m, as is appropriate for this example. 
We have also explicitly used the fact that = 1 in the coefficient of the 
mass term, and introduced the abbreviation 

da(rj) = S(l-J:rjl)Udrj k . 

We also introduce the notation for an average over a given by ((•)) = 
J(-)da(r))/ J da(r]). First we observe that (rjl) = (l/N) because every axis is 
treated identically. For all n > 2, we next observe that < ((rjk* — i]k) 2 ) — 
(l/N) as well, as follows from the previous equality plus the Schwarz in- 
equality. The implication of this fact is that there may be relatively little 
correlation between nearest neighbor values. [On the other hand, for n — 1, 
it follows that {(rjk* — Vk) 2 ) — (°>/N), indicating some significant correlation 
between nearest neighbor i] values.] We assume, and will later confirm, that 
the gradient term is generally the dominant term in the exponent. With that 
assumption it follows that the entire gradient term is of order N, i.e., 

i/€ 2 £( % * - Vk ) 2 a n - 2 ~ N . 

This relation holds for the simple reason that each of the N terms in that 
sum will be of order one just as would be the integrand for any single (TV- 
independent) integral. The value of k is well approximated by the result of 
a steepest descent approximation, namely as given by (N/ k) = na n ~ 2 where 
we have already used the (upper limit of the) estimate for (i] k * —i] k ) 2 and set 
N — 1 ~ N. We conclude that k ~ \ / Na 2 ~ n . With this choice for k it follows 
that the entire gradient term is 0(N) as desired. With these expressions we 
also learn that 

\m 2 K 2 a n ~ m 2 (Na 2 - n )a n ~ No 2 

which is of lower order than the gradient term justifying our neglect of it in 
performing the steepest descent approximation. For comparison of relative 
terms we may also examine a 'unit volume' where Na n = 1. In that case, 
when n > 3, the gradient term is 0(a~ n ) while the mass term is 0(a 2 ~ n ); 
for n = 2 the latter term should be interpreted as a = — In (ma) for some 
mass m. Thus for all n > 2 the mass term is divergent as iV — > oo although 
it diverges at a slower rate than the gradient term. Let us see what is the 
consequence of the divergence of the mass term. 

Suppose we want to consider the relation of S(h) for one mass term 
compared to that for another by means of a perturbation in the difference 
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Am 2 of the two masses. In particular, consider 

S(h) =C N j exp{nJ:h k ri k a n - |/t 2 S(%» - %) 2 a n " 2 

-\{m 2 + Am 2 )K 2 a"} k^" 1 dn da(r)) 
= C N J exp{Kj:h k r] k a n - |K 2 £(r? fc . - %) 2 a n ~ 2 

-\m 2 K 2 a n }(l - \Am 2 K 2 a n + • • •) /t^" 1 dn da(r]) . 

The point to observe in this calculation is that the term proportional to Am 2 
is, as compared to the first term of the series, proportional to Na 2 due to the 
presence of K 2 a n in the expansion. Thus we see that a perturbation in the 
mass is automatically divergent as N — > oo whenever n > 2, a fact already 
well known in conventional treatments of relativistic perturbation theory. 
Observe further that this divergence arises from the single variable k and is 
not a consequence of integrations over the many r\ k variables. 

Let us next consider the case of a quartic perturbation term, namely 

S(h) =C N J exp{Kj:h k r] k a n - \K 2 ^(r] k » - r] k ) 2 a n - 2 

— \mlK 2 a n — g n A 'Er) k a n } k^^ 1 dndairf) . 

Once again we assume to begin with that the gradient term is the most dom- 
inant term in the exponent, an assumption that leads to the same estimate 
for k, namely k ~ V Na 2 ~ n . The mass term is the same as before, while the 
new term is the quartic coupling which becomes 

g^4a n ~ g (Na 2 - n ) 2 NN- 2 a n ~ Na 4 ~ n . 

In this expression the factor 1/iV 2 in the middle arises from the estimate that 
(r] k ) ~ 1/N 2 . In fact, this value can be determined from 

, 2 2\ _ _ ^ + 2£h 
{VkVi ) - N ^ N + 2 ^ 

which holds exactly. Observe for n < 4 that the quartic term is strictly 
smaller than the gradient term, while for n = 4 the two terms are of the 
same order of magnitude. On the other hand, for n > 4 the quartic term is 
strictly larger than the gradient (or the mass term) and thus our assumption 
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of determining the behavior of k is unjustified. These three cases corre- 
spond exactly to the three cases of superrenormalizable, renormalizable, and 
nonrenormalizable fields in the usual perturbation analysis Here we see 
how these three different situations appear in the hyper-extreme spherical 
coordinates of the present paper. 

As an example of estimation obtained in this picture, let us determine the 
contribution to the vacuum polarization of the two-vertex bubble graph for 
n > 3 when estimated in hyper-extreme spherical coordinates. In particular, 
let us estimate the expression 

C N J(Z: n 4 r]t : a n ) 2 exp{-i/t 2 S(^ - %)V" 2 - \m 2 K 2 a n } k n ~ 1 dn da(r)) . 
Recall the usual prescription for normal ordering given by 

In terms of hyper-extreme spherical coordinates we see that 

4 4 4 2 a 2 2/ 2 2\ , ol 2 2\2 

: « Vk ■ = k Vk ~ 6« Vk( K Vk) + 3(« Vk) ■ 

Now we can use the fact that (ft 2 ?? 2 ,) — iVa 2 ""^ -1 = a 2 ~ n . All these terms 
then become of comparable value, and to the accuracy we are working we 
cannot distinguish a cancellation of a divergence (as with (: 4>l4>f :) = 0) 
from a noncancellation (as with (: </>| :: (f)f :) ^ 0). Rather we use the normal 
ordering to establish the proper combinatorics of the correlations. Thus 
we determine, at least to the correct order of magnitude, that the leading 
divergence of the expression given above, relative to the unperturbed term, 
is given by 

(Na 2 - n ) 4 a 2n £(^f) ~ N*a 8 - 2n a- n NN-* = Na 8 ~ 3n . 

In the middle estimate, and apart from k s , we have used the fact that the 
first sum would lead to a convergent integral were there a volume element 
(a _ra ), the second sum is translationally invariant (N), and the expectation 
of the 77 variables involves 1/iV 4 . The net result is proportional to the volume 
(N) as it should be. The remaining divergence is ultraviolet as it should be. 
It is readily confirmed that the result is appropriate for n > 3. 

The virtue of using hyper-extreme spherical coordinates is not necessarily 
for calculational purposes, but rather to gain an overview of where, in some 
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sense, the divergences that appear in quantum field theory originate. With 
the choice of hyper-extreme spherical coordinates we have exposed that the 
radius k — and especially its appearance in the integrand in the form k, n_1 — is 
the source of all the divergences. Let us make this point even more convincing 
by discussing what happens if that large power is removed! 

3 Divergence-free Quantum Field Theory 

For the sake of discussion let us replace the term that we have identified as 
the source of the divergences by something harmless. In particular, let us 
replace k n ~ 1 by 1, i.e., effectively, in this term and in this term alone, setting 
N — 1. This is a very drastic change and it has a number of consequences 
not the least of which is the loss of a local interaction with which to make 
a relativistic model. We will comment on this very point later, where we 
will briefly recall the principal argument, which is presented elsewhere ||, 
in which an analogous replacement appears in such a way that locality in 
the continuum limit is retained. For the present — and this is the main point 
of the present paper — we wish to emphasize the profound consequences that 
arise when such a change is made in a lattice space functional integral. 

With the indicated change in the integration measure, the lattice space 
generating functional of interest becomes 

S'(h) = C' N J e W {KY,h k n k a n - i K 2 7S( % . - n k fa n ~ 2 - ±mjfc 2 E^ a " 

-g Q K p J:r] p k a n } dn 5(1 - S^) n dn k . 

One immediately sees that the integration over the single variable k is normal 
and without any special feature. On the other hand, the parameters of the 
model, i.e., Y, mo, and go, need not behave as they previously did, and 
indeed, the cutoff dependence of these parameters is one of the interesting 
features that needs to be determined. To that end let us first examine the 
special case in which Y = = go and only m remains. In that case 

S'(h) = C' N J exp{Ki:h k r] k a n - lm 2 K 2 a n } dn 5(1 - T,rg) II dn k . 

Although « > 0, the symmetry of the measure under n — > —rj lets us treat 
the integral over k as a simple Gaussian integral, and the result (changing 
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" N as needed) becomes 



S'(h) =C' N J exp{ Imfim^) 2 ^ } 5(1 - E V l) Udrj k 



The remaining integrals over the 77 variables may be approximated by simply 
replacing the factor rj k rji by its average value of 6m/ N. When this is done the 
result becomes 



and in the continuum limit the only way for the result to be nontrivial is that 
tUq = m 2 /N, where m 2 is a cutoff independent quantity. This dependence of 
the bare mass is of course quite different than in the usual case. Furthermore, 
it is noteworthy that a perturbation in the mass term leads to a series that 
is order-by-order finite. In particular, consider the expression 



S'(h) =C' N J e^inT.hk^ - \{m 2 + ^m 2 )N- l K 2 a n } dn6{\ - T 1 r 1 2 )Yidr lk 

= C' N J exp{Ki:h k r] k a n - \m 2 N~ x K 2 a n } 

x (1 - \/\m 2 N~ l K 2 a n + ■ ■ •) dK 5(1 - Y.rg) II dr] k . 



The essential fact to observe about this equation is that the term of order 
Am 2 is actually 0(1) due to the fact that k 2 is large of order Na~ n . Thus 
there is no divergence in the perturbation series as long as we choose the 
mass renormalization as indicated. 

Next let us introduce a nonlinear coupling term. This leads to the ex- 
pression 



In view of the behavior of the integral over n, it follows that this expression 
will be well defined and contributing if and only if 



S'(h) = exp{ lm 2 N- 1 m 2 a n } , 





S'(h) = C' N J exp^iV^E/^a™/ 2 - I m V 

-g (fN^-^ydQ 6(1-^1) II d Vk , 
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where, as usual, the constant C' N has been readjusted to account for the 
change of variable. A perturbation expansion in the nonlinear coupling con- 
stant leads to 

S\h) = C' N J exp{gN 1/2 J:h k ri k a n/2 - \m 2 q 2 } 

x (1 - gg^N^Eril + ...)dg5(l- S % 2 ) II d Vk , 

It follows, relative to the first term, that the first-order correction is 0(1). 
A similar analysis shows that each of the terms in the perturbation series 
is finite. In other words, the perturbation series is order-by- order finite. 
It is particularly noteworthy in this regard that go is small (specifically 
go oc N^a^l 2 ^™ ), meaning that the factors that it multiplies need to be 
allowed to become large in order for the nonlinear term to contribute. A large 
value for the factors multiplying the coupling constant is analogous to what 
happens in the usual case where the factor n N ~ l is present in the integrand. 
But not only are these large, they also have a broad range of values, and 
that situation is quite unlike what happens in the usual case when the factor 
k n ^ is present. 

Finally we turn our attention to the gradient term. For this purpose we 
first let m = = go so that we may focus on the essentials. The expression 
of interest reads 

S'(h) =C' N J exp{ K m kVk a n - |/€ 2 r£( % * - Vk) 2 a n - 2 } dn8{\ - £ % 2 ) Ud Vk . 

Clearly the role of the gradient term is to introduce correlations into the direc- 
tion field that have been, in the absence of that term, lacking. Nevertheless, 
since rjl ~ 1/N, we observe, as previously noted, that < (r] k * —r] k ) 2 ~ 1/N 
as a consequence of the Schwarz inequality. Thus to make this term a con- 
tributor it is necessary that 

K 2 YnVk*-rik?a n - 2 ~l 

a condition that requires that 

Y ~ {Na- n )- 1 a 2 - n {T 1 {r ]k *-r ]k ) 2 )- 1 ~ N^a 2 . 

Finally, in this regard, we combine the several terms needed to make up 
the full theory (minus, of course, the factor k n ~ 1 in the integrand!). Based 
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on the fact that each of these terms has a bounded contribution, it follows 
that we may choose as our basic lattice-space expression 

S'(h) = C' N J exp{/t£/i fe% a n - \K 2 N- l a 2 ^ k * - i] k ) 2 a n - 2 - \m 2 N- l K 2 a n 

-gN- 1 a W2 - 1)n K p Tlr] p k a n } dn 5(1 - E^) n dr] k . 

This equation, with the indicated cutoff-dependent coefficients, represents 
our proposal for obtaining a nontrivial (and generally non-Gaussian) con- 
tinuum limit in the special case that the factor /t^ -1 is removed from the 
integrand. 

Moreover — and this is of central significance — the form of the renormal- 
ization of coefficients involved is entirely multiplicative and in no way in- 
volves normal ordering. This has the general consequence that for a covariant 
quantum field theory defined in some way which is analogous to the present 
discussion, subtractive renormalizations would not be required; instead, all 
required renormalizations would be strictly multiplicative. 

Let us now turn our attention to a very brief discussion of how we believe 
a scenario that is schematically similar to that of the present paper may well 
be adapted to a relativistic self-interacting scalar field. 

4 Application to Relativistic Models 

The quantization of any system entails ambiguities including the introduction 
of additional potential terms to the Hamiltonian operator or to the action 
that are proportional to h. In field quantization it is generally acknowledged 
that some such additions are necessary. The proposal we have in mind adds 
an auxiliary, /i-dependent potential P to the lattice action designed to dra- 
matically alter the continuum limit. Specifically, the potential P is designed 
to make the ground-state distribution of the field at a sharp time have a gen- 
eralized Poisson distribution in the continuum limit. One way to arrange this 
is for the distribution to be chosen already on the lattice to have a generalized 
Poisson distribution. To illustrate this proposal we choose the characteristic 
function for the lattice-space, sharp-time field distribution 1^(0) | 2 to have 
the form given by 

J 9k<t>ka s | \[r(0) | 2 n' d(f) k = exp{-/[l - e lEWfcaS ]p(0) n'd<M 
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where p > 0. Here, gk denotes a test sequence at a constant lattice time, 
and the sums and products run over the set of lattice sites at a constant 
lattice time. To ensure that the model has a unique translationally invariant 
state it is necessary and sufficient that / p(0) H' d(j)k = oo, which implies that 
we deal with a generalized Poisson distribution |4[]. This divergence can be 
traced to a necessary singularity in the integrand for small field values which 
is accommodated by the ansatz that 



n' fc [s;/3 fc _^ 2 ]7 ■ 

Here J 2 is a normalization constant, W(4>), with W(0) = 0, accounts for 
the intermediate and large field behavior and leads to the finiteness of the 
moments, {/3k}, k a lattice point, is a set of positive constants (e.g., /3& oc 
exp[— E| =1 |/cj|] ) subject to the normalization condition = 1 where the 
sum and product are over all N' = (2L + l) s lattice sites at a constant 
lattice time and with periodic boundary conditions, and finally 7 = 1/2 + 
£/(2iV'), < £ < 1. With this implicit ansatz for the lattice ground state 
the lattice Hamiltonian is defined by 



d 2 1 d 2 ^(4>) 



Although it is generally difficult to determine an analytic expression for the 
ground state ^(0), it is possible to argue that in the continuum limit the local 
form of the auxiliary potential P is proportional to <f)(x)~ 2 with a coefficient 
proportional to fi 2 . By using the freedom to choose W we hope to ensure 
that the remainder of the potential V has the necessary lattice form for a 
covariant continuum limit. Of course, once the proper Hamiltonian is in 
hand, it is possible, in principle, to build various correlation functions. 

We finally comment on the relation between the proposal of the present 
section for a relativistic field and the discussion in the foregoing sections. 
That connection is readily made when we once again introduce hyper-extreme 
spherical coordinates, this time on a spatial slice of the lattice. In particular, 
the form for the characteristic functional for the sharp-time field becomes 

exp{-2J 2 /[l - e «s's«a'] e -2Wh) n'jE^^^ 2 ]-^- 1 ^^ 5(1 - S^ 2 )n'rf % } . 
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Observe that the factor k 1 ^ is a result of a quotient between a factor k n ' 1 
in the numerator which arises from the measure and a factor 
that appears in the denominator. The cancellation of the large power k n that 
arises is exactly analogous to the dropping of the factor k n ^ in the previous 
section of this paper! Lastly, we note that the choice of the resultant power, 
i.e., —1 — £, is designed to give a divergent result for / p(<p)Il'd<f) k , and yet, 
with the help of W, make every moment of the smeared, sharp-time field 
T,'g k (f) k a s = KZ'g k r] k a s finite. 

Dedication 

It is a pleasure as well as an honor to dedicate the present article to the 
memory of Hiroomi Umezawa, whose warmth and enthusiasm in all matters, 
scientific and otherwise, was universally recognized and appreciated. 
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